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Abstract 

Lie-Poisson structure of the Lorenz'63 system gives a physical in- 
sight on its dynamical and statistical behavior considering the evo- 
lution of the associated Casimir functions. We study the invariant 
density and other recurrence features of a Markov expanding Lorenz- 
like map of the interval arising in the analysis of the predictability of 
the extreme values reached by particular physical observables evolving 
in time under the Lorenz'63 dynamics with the classical set of param- 
eters. Moreover, we prove the statistical stability of such an invariant 
measure. This will allow us to further characterize the SRB measure 
of the system. 
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1 Introduction 

In 1963 E. Lorenz by a drastic truncation of fluid-dynamics equations govern- 
ing the atmospheric motion obtained a system of ODE which he proposed as 
a crude yet non trivial model of thermal convection of the atmosphere [Lj. As 
a matter of fact, the Lorenz model is today understood as a basic toy-model 
for the evolution of Earth atmosphere's regimes, like zonal or blocked circula- 
tion or climate regimes (e.g. warm and cold), in which dynamics is described 
by equilibrium states [CMPj . [Se] . In his work Lorenz showed such system 
to exhibit, for a large set of parameters values, a peculiar chaotic behavior, 
that is exponential sensitivity to perturbations of initial conditions and the 
existence of a global attracting set for the flow nowadays called generalized 
nontrivial hyperbolic attractor. Although there exists an extensive literature 



on the subject, we refer the reader to Sp for a rather comprehensive overview 



on this problem and to [V] for a recent account on the progress made on the 
rigorous analysis of the Lorenz '63 ODE system and the relationship between 
this and its more abstract counterpart, the geometric Lorenz model, intro- 
duced in the second half of the seventies ([G]. |ABS] and |GW[ ) to describe 
the geometrical features a dynamical system should posses in order to exhibit 
the same asymptotic behavior as the Lorenz one. An affirmative answer to 
the long-standing question whether the original Lorenz '63 flow fits or not 
the description of the one modeled by the geometric Lorenz, which means it 
supports a robust singular hyperbolic (strange) attractor, has been given by 
W. Tucker in [T] by means of a computer assisted proof. As a byproduct, 
Tucker also proved that Lorenz flow admits a unique SRB measure supported 
on the strange attractor. Further results about the characterization of the 
set of geometric Lorenz-like maps are given in |LM[ . 

In 2000 it has been emphasized that Lorenz '63 model and the Kol- 
mogorov one, considered as a low-order approximation of the Navier-Stokes 
equations, belong to a particular class of dynamical systems, named Kolmogorov- 
Lorenz systems |PPlj , whose vector field admits a representation as a sum of 
a Hamiltonian SO (3)-invariant field, a dissipative linear field and constant 
forcing field (see also |PP2] for an extension of this analysis to the Lorenz 
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'84 model). Moreover, they proved that the chaotic behavior of these models 
relies on the interplay between dissipation and forcing. 

More specifically, and more recently, in |PM] it has been shown that 
the effect of the dissipative and forcing terms appearing in the previously 
described decomposition of the Lorenz '63 vector field, with the classical set 
of parameters, is to induce chaotic oscillations in the time evolution of the 
first integrals of the Hamiltonian system associated to the Lorenz '63 model, 
namely the Hamiltonian and the Casimir function for the (+) Lie-Poisson 
brackets associated the so (3) algebra [MRj . which represents the angular 
momentum of a free rigid body in the Kolmogorov-Lorenz representation of 
geofluid dynamics introduced in [PP.1]. In particular, it has been shown that 
two subsequent oscillation peaks in the plot of the Casimir function C as a 
function of time are related by a map $ of the interval similar to the one 
originally computed by Lorenz in [L] depicting the functional dependence 
of two subsequent maximum values assumed by the third coordinate of the 
flow as a function of time. We remark that, being C the square norm of the 
flow, the similarity of the plots of this two maps is therefore not surprising. 
In [PMj, the recurrence properties of $ are also studied, which allows to 
characterize the trajectories of the system through the number of revolution 
they perform around the unstable point lying on one side of the plane x + y = 
when the initial condition is chosen on the opposite side ( |PM| figg. 2,10 
and 11). 

In our paper we clarify what stated in |PM] by giving an account in 
the first section of the rigid body formulation of the Lorenz '63 model and 
constructing, in the next section, a Markov expanding Lorenz-like map T 
of the interval being the reduction to [0,1] of $. Both maps are in fact 
derived throughout the Poincare map associated to the surface in the config- 
uration space of the system corresponding to the set of maxima reached by 
the Casimir function during its time evolution. Hence, we will study the in- 
variant measure under the dynamics defined by T characterizing its density 
and consequently the SRB measure of the system. Furthermore, we anal- 
yse the recurrence properties of the dynamics induced by T clarifying more 
rigorously what stated in Section IV B of |PM] . 

We will also perturb the system by adding an extra forcing term which 
will eventually cause the system to loose its symmetry under the involution 
R : (x,y, z) —> (—x,—y,z) of 1R 3 . Due to the robustness of the attractor, 
i.e. persistence under perturbations of the parameters, proved in [Tj, maps 
analogous to T can be defined and studied, and their statistical properties 
analysed as in the unperurbed case. Therefore, such perturbation of the 
Lorenz '63 field will only have the effect to induce a change in the statistics 
of the invariant measure for the system, the SRB measure. We will prove 
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that this change could be detected by looking at the deviation of the in- 
variant density of the perturbed map with respect to the unperturbed one. 
Such result would confirm what has been empirically shown in [CMP] about 
the impact of anthropogenic forcing to climate dynamics of the northern 
hemisphere. 

We also believe that this analysis could also be pursued in the case of 
more general iV-dimensional models such as those introduced by Zeitlin in 
[Z] to approximate in the limit of N tending to infinity the dynamics of the 
atmosphere in absence of dissipation and forcing. 

We will present elsewhere our contributions in these directions; here we 
prove the first non-trivial result about the statistical stability of the invariant 
measure for T. The technique we propose is new and we believe it could be 
applied as well for other maps with some sort of criticalities. 

We remark that, in particular, the distribution of the return times of a 
measurable subset of [0, 1], which can be derived directly from the invariant 
measure of T, could be useful in studying the statistics of extreme meteoro- 
logical events. 

2 Return Lorenz-like maps 

2.1 Rigid body formulation of the Lorenz '63 model 

It can be shown ( |PPlj ) that the Lorenz '63 ODE system [L], 

Xi = —axi + ax 2 

x 2 = —xix 3 + pxi - x 2 (1) 
X 3 = X\X 2 - /3x 3 

can be mapped, through the change of variables 

U\ = X\ 

u 2 =x 2 , (2) 

U3 = x 3 - (p + cr) 

to the ODE system 

Ui = —au\ + 0U2 

u 2 = —U1U3 - au\ - u 2 (3) 
ii 3 = uiu 2 - (3u 3 - (3 (p + a) 

representing the evolution of a Hamiltonian system whose configuration space 
is the 5*0 (3) group, subject to dissipation and to a constant forcing. That 
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is, denoting by 

{F, G} := col ad^ G x) = x ■ VF x VG (4) 

the Lie-Poisson brackets associated to the symplectic 2-form u\ defined on 
the cotangent bundle of SO (3) [MR], @ reads 

v-i — {ui, H} — (Au) i + fi , i = l,2,3, (5) 

where: 



H (u) := -u ■ Qu + h ■ u (6) 

is the Hamiltonian of a rigid body whose kinetical term is given by the 
matrix Q := diag (2, 1, 1) , while h := (0, 0, —a) is an axial torque; 

• A := diag (a, 1, 0) is the dissipation matrix; 

• / := (0, 0, —(3 (p + cr)) is a forcing term. 

This representation allows to study the Lorenz system as a perturbation 
of the Hamiltonian system 

Vi(u):={ui,H}, 1 = 1,2,3, (7) 

admitting, as in the case of a rigid body with a fixed point, two independent 
first integrals H and the Casimir function C for the Poisson brackets Q 
|MRj . In fact, when rewritten in this form, it follows straightforwardly that 
the system is non chaotic for a = |PP1] while, for a ^ 0, the values of C 
and H undergo chaotic oscillations |PM] . 

Moreover, when passing to the representation ([5]), the symmetries of the 
system are preserved as well as other features such as the invariance of the 
£3 (^3) axis and the direction of rotation of the trajectories about this axis. 
The critical points of the velocity field of the system are then 



ci := 



VJIp- 1)> VHp~ 1)> -{?+ 1)) , c 2 := {-xi (ci) , -x 2 (ci) , x 3 (ci)) 



and c := (0, 0, - (p + cr)) . 

We also remark that (|5| can be rewritten in the form 



u = v — w , (9) 
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where v is the divergence free field ([7]) and 

R 3 9m — >w(u) := Au-f = S7K(u) Gi 3 (10) 

with 

K{u) := -u- Ku- f -u (11) 

a convex function on IR 3 . Notice that the fields v and w are orthogonal in 
L 2 (r£>;IR 3 ) , where B := {u G R 3 : ||u|| < 1} is the unitary ball in R 3 and 
rB denotes the ball of radius r. 

The decomposition of the velocity field as the sum of a divergence free 
field ad a gradient one, together with the appearance in the Hamiltonian 
description of the flow of the Lie-Poisson brackets Q in the space reference 
frame of the rigid body, i.e. right translation on SO (3) , is standard in fluid 
dynamics [MRJ and can be seen as another source of analogy between 
the Lorenz '63 model and Navier-Stokes equations |PPlj . |FJKTVj . 

2.2 The return map on the set of maxima of the Casimir 
function 

/II 

let 



If Uq is any non stationary point for the field (3) such that ||w || < T/^p with 



min{t e speck}, and C (t) := w )| 
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m:=MC{t) ; M:=supC(t) . (12) 



Clearly, m > since C > 0. Moreover, M < oo since it has been shown in 
jPPT] that C (t) < ^ where, for our choice of parameters, 

= /3(p + a) . (13) 



VAa 

To construct the function which links two subsequent relative maximum 
values of C (t) we proceed as follows: 

• first we identify the manifold S in the configuration space of the system 
corresponding to the relative maxima of C (t) , 

• then we construct a map of the interval [0, 1] to itself as a function of 
the map of the interval [m, M\ of the possible values of C (t) in itself, 
which can be defined through the Poincare map of this manifold. 

The existence of the aforementioned Poincare map follows from the exis- 
tence of the return map computed by Tucker in [T] . 

Throughout the paper we will assume a — 10, p = 28, /3 = |. 
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2.2.1 Identification of S 

By Q we get 



C (u) = -2 



E(u) 



(3{p + a) 



with 



Therefore, 



E (u) := au\ + u\ + (5 ( u 3 



(p + cr) 



£ : = <! u G M 3 : C (u) = S- = <J u G M 3 : E (u 



0(p + a) 



(14) 



(15) 



(16) 



as already noticed in |PMj . is an ellipsoid intersecting the vertical axis (u 3 ) 
in the origin and in cq. This also implies, M = p + a. Clearly, c\, c 2 G £. 
Moreover, by ([14} and Q 



(7 (u) = 2V-E -[«xVif + ViT] (u) 



(17) 



a (a - 1) + (0 - 1) « 3 + 



p + cr \ p + a 



+/3 2 «3 + 



p + cr 



0- 



,p + a 



u 3 + 



p + cr 



Let us set z := w 3 + then 



116K 3 : a 2 ^ + u 2 , - 



£' := \ u G M 3 : C (u) = 



a (a - 1) + (/3 - 1) z + 



p + cr 



U\U 2 



+P 2 z z + 



p + cr 







We remark that, denoting by R the involution 

1R 3 3 u = {u\,u 2 , Us) 1 — > Ru := (—Mi, — u 2 , u 3 ) G 1R 3 

leaving invariant the field it, RE = £ and R£' = £'. 
Consider the diffeomorphism 



(19) 



qi = O (z)ui, i = l,2 ; q 3 = z 



(20) 



such that for any fixed value of z, (z) is an orthogonal matrix diagonalizing 
the symmetric quadratic form 



a(a-l) + (p-l)z + 



p + o 



C1C2 , (Gt 2 , 
(21) 



namely, setting A (z) = O l (z) diag (Ai (z) , A 2 (z)) O (z) , 

u ■ A (z) u = q ■ O (z) A (z) O t (z) g = Ai (z) q\ + A 2 (z) q\ . (22) 

with 



Ai(*) 



a 2 + l + J (a 2 -l) 2 +[^ + a{a-l) + {(3- 1) z 



a 2 + l- J (a 2 - i f + [ep + <? (o- - 1) + W - 1) z 



Under this change of variables 



£' = {9 e : Ai (g 3 ) 9i + A 2 (q 3 ) q\ + /3 2 g 3 (93 + 



p + a 



(23) 
(24) 

(25) 



Since Ai (q 3 ) is positive for any choice of the parameters (3,p,a and q 3 , the 
equation giving the intersection of £' with the planes parallel to q 3 = 
(w 3 = — ) can have a solution only if A 2 (#3) is negative, that is for 



93 > 



93 < 



^-3) + ^ 

Cr(cr+1) + ^ 
5-1 



W 3 > 



U3 < 



1 



- 1 



a (a - 3) + p 
<t(<t + 1)+P 



(P+g) 
2 

(p + a) 



; (26) 

• (27) 



Therefore, for g 3 7^ (w 3 7^ — ^jr^), these intersections are hyperbolas while, 
if g 3 = or g 3 = — (it 3 = — (p + cr)), from the definition of £' we get 
the equations 



if 93 = 0, 



2 2 1 2 

(T ti 1 + U 2 — 



(p + a) 



+ a(a-l) 



u x u 2 = ; 



(28) 



• if 93 = 



(p±g) 
2 ' 



2 2 1 2 
cr w x + u 2 — 



a(a-l)-{(5-2) 



P + cr 



U\U 2 = . 



(29) 
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Since for our choice of the values of the parameters of the model, 



A 2 (0) = 



+ a(a-l) 



< 



A 2 



(p + a) 



a 2 + 1 _ u a 2_ 1 y + < T (< T _i)_(0_2) 



(30) 

< 
(31) 



the intersection of 5' with the planes z = q% = and z = g 3 = ^±£0 are the 
straight lines. The manifold in i? 3 corresponding to the relative maxima of 
C (t) is 



E := \u G 



{« G M 3 : C(u) = , C7(«) < o} 



(32) 



p3 . 



l(P+^) 2 
4 



: < " M ' : < a 2 u\ + ul-[ a (a- 1) + (f3 - 1) (u 3 + + «i«2+ 
+/3 2 (% + ^) («3 + P + ^)<0 

Since for our choice of the parameters 



a (a - 3) + p 



(p + a) 



> p + a , 



(33) 



E is composed by two closed surfaces in R 3 , E + and E_, such that i?E + = E_ 
and intersecting only in the critical point cq. 

By definition, Vu G £, the vector -u (w) is orthogonal to the vector VC (u) , 
hence it belongs to the plane spanned by VE (u) — (VE ■ VC) (u) VC (u) and 
(VC x V£) (u) , where 



u 2 [(I3-I)u 3 + P<¥] 
VCxV£ = 4[ - Ul [(/3-a)u 3 + p^] 

— (a — 1) 



(34) 



and, since C is a constant of motion for the Hamiltonian field v, 
Vw e£, (w- VC) (u) = 0. 
Moreover: 



• \VE\ ]s, |VC| ]s and |VC x V-E| ]s are always different from zero; 
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from (Tub and d32b 



C if) = (u ■ WCj (t) = [ u ■ V 



-2 



(p + a) 



(*) , (35) 



hence, Vw G c?E, u (u) is parallel to VC x V-E, that is tangent to E. 
Therefore, u is transverse to £\<9£ and since Vm G E\SE, 
(7 (u) = (u ■ VC") (u) = -2 (u ■ WE) < (u • V£) > , (36) 
the direction of -it (u) points outward the bounded subset of IR 3 , 

,(P + ^) 2 



it G IR 3 : £ (u) < /3- 



(37) 



2.2.2 Parametrization of E 

If r G (0, p + a) , 7 := rB fl £ is a regular closed curve. Therefore, we can 
parametrize E + choosing an appropriate arc of 7 as coordinate curve of the 
parametrization, that is there exist an open regular subset Q of IR 2 and a map 
b + G C 1 (n, IR 3 ) n C (H, IR 3 ) such that 



ui = bf (y, z) 

u 2 = 63 (2/, z) , (y,2)6fl. 



(38) 



Moreover, we can choose the parametrization such that z = r 2 , therefore the 
coordinate curves b\ (y) satisfy the equations 



C(bt (y)) = z 

c(K (y)) = o 



(39) 



We remark that the tangent field to 7 is parallel to VC x WE, while, if 
C denotes the coordinate curve 6+ (z) , the tangent field to ( is parallel to 
WE x (VC x WE) . 

Similar arguments also hold for £_. Furthermore, 



Q3(y,z)^b (y, z) := Rb + (y, z) G 



(40) 



is easily seen to be a parametrization of £_ sharing the same properties of 
6+ ' 



10 



2.2.3 Return maps on £ 

The evolution of the system maps in itself the ball (3 (p + a) B, £ C /3 (p + a) B 
and the velocity field is transverse to E\<9£ and £\£. 

For our choice of parameters p, f3 and a, it has been shown in [T] that there 
exist periodic orbits crossing a two-dimensional compact domain A contained 
in the plane n := {u G 1R 3 : u 3 = 1 — (p + cr)} , which is also intersected by 
the stable manifold of the system W£ along some curve r . Furthermore, the 
first eight shortest periodic orbits have been rigorously found in |GTj . Notice 
that by symmetry if p(t,u) , u G £+, is a periodic orbit, Rtp (t, u) is also a 
periodic orbit and either Rip (t, u) = (p (t, u) or Rp (t, u) = p (t, Ru) , where 
Ru G that is periodic orbits are either symmetric or appear in couples 
whose elements are mapped one into another by R, as already remarked in 

Since A is easily seen to be contained in 

|w G R 3 : C[u) < o} nvr , (41) 

these periodic orbits then, necessarily cross £ which is also possibly inter- 
sected by along some curve V lying in the half-space 

{u G M 3 : u 3 > 1 - (p + cr)} . (42) 

Therefore, if Uq G S\r lies on a periodic orbit of period t , there exists 
an open neighborhood N 3 uq and a C 1 (N, M) map r such that r (uq) = to 
and p T ( u ) (u) G E for any u G N. Then, 

NO J:\T3ui — > P E (u) := p r[u) (u) G E . (43) 

Moreover, it has been proved in [T] that A\r is forward invariant under the 
return map on 7r and that on A there estists a forward invariant unstable 
cone field. These properties are also shared by a compact subset A' C E such 
that any open subset of A' is diffeomorphic to a open subset of A. Hence, 
Ps admits an invariant stable foliation with C 1+L , i G (0, 1) leaves. 

2.2.4 Construction of the map T 

Let pW := Pe ± . By the parametrization previously introduced for E + , there 
exists an open subset Q' C Q\T', with V := (b+y 1 (T) , and a C 1 (Q',R 2 ) 
map 

Q'3{y,z)^S(y,z)en' . (44) 

such that, V (y, z) G ft', 

(b+oS)(y,z):=(P^ob + )(y,z) , (45) 
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Furthermore, 

G(y,z):= (Cob + oS)(y,z)=0. (46) 

Let Si, S 2 be respectively the first and the second component of S. Since b + 
is a diffeomorphism and the components of VP = VC are different from 
zero on E + , V(y,z) G ft', d y G(y,z) ^ 0. Thus, by the implicit function 
theorem, V(yo,z ) G ft', there exist two open interval (2/1,2/2) , (zi,z 2 ) such 
that (y ,z Q ) G (2/1,2/2) x (zi,z 2 ) Q ^' and a unique C 1 ((z l} z 2 )) map 

(zi,z 2 )3 z^y.= U{z)e( yi ,y 2 ) (47) 

such that, G (zi,z 2 ), G (U (z),z) = and, V(y,z) G (2/1,2/2) x (zi, z 2 ) 
such that y^U(z), G (y, z) ^ 0. 
Therefore, let 

(*!, z 2 ) 9 z ^ F (z) := S 2 (U (z) , z) G (zi, z 2 ) . (48) 

Notice that, since b + G C 1 (ft) , 5 = (6+) _1 o P« o 6+ is C 1 (ft') if and 
only if p( + ) is, and so are U and V. 
Moreover, by symmetry, 

b- oS = Rb + oS = PP (+) ob + = PP (+) o Rb~ = P { -) o b~ . (49) 

Hence 5 = (ft") -1 o pH 

Clearly, [zi,z 2 ] C [m V (r*) 2 , p + tr] , with r* := inf{r > : rPnS ^ 0}. 

Let m + G S + ,(y + ,2; + ) G ft be such that P^ (u + ) = c and 
b + (y + , z + ) = u + . Setting 

[zi, z 2 \ 3z ^X{z):= G [0, 1] , (50) 

z 2 - Zi 

we define 

[0, 1] 3 s 1 — > T (s) :— X oV o X' 1 (s) G [0, 1] . (51) 

Hence, by construction, T is a C 1 ((0, 1) \ {xq}) map, where xo := X (z + ) , 
and, since there exits 1 G (0, 1) such that Ps admits an invariant stable 
foliation with C l+i leaves, then T is also C l+i ((0, 1) \ {x }) . 

3 The invariant density for the evolution un- 
der T 

In this section we compute the density of the unique (by ergodicity) abso- 
lutely continuous invariant measure for the map T and we prove its statistical 
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stability. For the construction of the density we use the techniques recently 
introduced in the paper |CHMVj (see also |BHj for results related to similar 
maps), which dealt with Lorenz maps admitting indifferent fixed points be- 
sides points with unbounded derivative. For the statistical stability we will 
follow the recent article |BV] . but with some new substantial improvements. 
The techniques used in |CHMVj turned around Young's towers |LSYj and, 
substantiated by a careful analysis of the distortion, led to a detailed study 
of the density of the absolutely continuous invariant measure, of the recur- 
rence properties of the dynamics and of limit theorems for Holder continuous 
observables. This analysis could in particular be carried over when the map 
has a derivative larger than one at the fixed point, as in the case we are going 
to treat, but possibly smaller than one in some other point (see below). We 
remind that whenever the Lorenz map is a Markov expanding map with fi- 
nite derivative, it could be investigated with the spectral techniques of Keller 
[K] . Young's towers are useful when the map looses the Markov property, 
but preserves points with unbounded derivative. This has been analysed in 
[KDO], see also |OHL] when there are critical points too. Our main effort 
will be in investigating the smoothness of the density. We will show that 
such a density is Lipschitz continuous on the whole unit interval but in one 
point. The argument we produce is a strong improvement with respect to 
the result achieved in jCHMVj (and applies to it as well), where we solely 
proved the Lipschitz continuity on countably many intervals partitioning the 
unit interval. We stress that, as far as we know, this is the first result where 
the smoothness of the density for Lorenz like maps is explicitly exhibited. 

Notations: With a n ~ b n we mean that there exists a constant C > 1 
such that C~ l b n < a n < Cb n for all n > 1; with a n ~ b n we mean that 
lim^oo = 1. We will also use the symbols "O" and "o" in the usual sense. 

The analysis we perform in this section applies to a large class of Lorenz- 
like maps which includes in particular those whose behavior is given by the 
theoretical arguments of the preceding section and by the numerical investi- 
gations of the paper |PM] . The map T (fig.l) has a left and a right convex 
branches around the point < Xq < 1; the left branch is monotonically in- 
creasing and uniformly expanding even at the fixed point 0, while the right 
one is monotonically decreasing with the derivative bounded from below by a 
constant less than one; at the cusp, located at xo, the left and right derivative 
blow up to infinity. Both branches are onto [0, 1] and this makes the map 
Markovian. Moreover, we recall that our map is C l on [0, l]\{xo} and C 1+t , 
i E (0,1), on (0,x ) U (x ,l). 

The local behaviors are (c will denote a positive constant which could 
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Figure 1: Normalized Lorenz cusp map for the Casimir maxima 



take different values from one formula to another): 

T(x) = a'x + (3'x 1+ ^ + o(x 1+ ^); x 0+ . . 

DT(x) = a' + cx^ + o(x^), a' > 1; > 0; V > 1 ' 

T(x) = a(l -x) + (3(1 - x) 1+K + o((l - x) 1+K ); x 1~ 
DT(x) = -a- c(l - x) K + o((l - :r) K ), < a < 1, /? > 0, k > 1 ' 

(53) 

= 1 - A'(z - a;) B ' + o((x - rc) B '); x -> Xq, A' > , , 

= c(x - x) 5 '- 1 + o((x - x)^- 1 ), < B' < 1 ' ( J 

T(x) = 1- A(x- x ) B + o{{x - x ) B ); rr ->■ x+, A > , . 

DT(i) = -c(x - xo)*- 1 + o((x - xq) 5 - 1 ), < B < 1 ' 

We set .B* := ma,x(B, B'); moreover we set T\ (resp. T 2 ) the restriction 
of T to [0, xq] (resp. to [rr ,l]). A key role is played by the preimages of 
x since they will give the sets where we will induce with the first return 
map; so we set: a := T 2 ~ 1 x ; a' := T 1 ~ 1 x ; a' p = T^ p a' ; a p = T 2 1 T 1 *- p 1 - ) a' , 
p > 1. We also define the sequences {b p } p >i C (x ,a ) and {O p >i C (oq,x ) 
as Tb' p = Tb p = a p _i. The idea is now to induce on some domain / and to 
replace the action of T on I with that of the first return map Tj into I. We 
will see that the systems (/, Tj) will admit an absolutely continuous invariant 
measure fij which is in particular equivalent to the Lebesgue measure with 
a density pi bounded from below and from above. There will be finally a 
link between the induced measure \ij and the absolutely continuous invariant 
measure \i on the interval, which will allows us to get some informations on 
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the density p of pi. The principal set where we will choose to induce is the 
open interval / = (a' , a )\{a;o}- The subsets Z p C I with first return time p 
will have the form 

Z\ = (a' ,&i) U (&i,a ) (56) 
Z P = (b' p _ 11 b' p )U(b P X-i) P>1- (57) 

We will also induce over the open sets (a' n , a' n _x) and (a n , a n+ i), n > 1, simply 
denoted in the following as the rectangles I n . In order to apply the techniques 
of |CHMVj . we have to show that the induced maps are aperiodic uniformly 
expanding Markov maps with bounded distortion on each set with prescribed 
return time. On the sets I n the first return map Tj n is Bernoulli, while the 
aperiodicity condition on / follows easily by the inspection of the graph of the 
first return map T T : I — > I showing that it maps: (a'^b^) onto (x ,a ); the 
intervals (^,5' /+1 ), I > 1, onto the interval (a' ,x ) and (&i,a ) onto (x ,a Q ). 
Finally, Tj sends the intervals (bi + \,bi), I > 1 onto (a' ,xo). Bounds on the 
distortion of the first return map on I and on the I n can be proved exactly in 
the same way as in Proposition 3 of |CHMVj (we defer to it for the details) 
provided we show that the first return maps are uniformly expanding The 
proof of this fact is given in the next Lemma and it requires a few assumptions 
which can be checked numerically with a finite number of steps and by a 
direct inspection of the graph of T. With abuse of language we will say that 
the derivative is larger than 1 if its absolute value is larger than 1. 



Lemma 1 Let us suppose that in addition to (52)- (55) the map T satisfies 
the assumptions: 

(i) d(i,o) : = inf xe ( fel ,ao) \DT(x)\ > 1; 

(ii) |OT(6!)| > OT(oJ,); 



1 The Lorenz-like map considered in CH MVj was C 2 outside the boundary points and 
the cusp; our map is instead C 1+ ' . This will not change the proof of the distortion in 
[CHMV] and all the statistical properties which follow from it. As a matter of fact, in the 

initial formula (5) in |CHM Vj . we have to replace the term ^tH) \T q ( x ) ~ T q (y)\ with 

— — 'T 9 (x) — T q (y)\ L , where Ch is the Holder constant larger than depending only on 



\DT((,)\ 

T and £ is a point between the iterates T q (x), T q (y). The only delicate point where the 
C 1+c assumption could give problems is the summability of the series at point (i) in the 
statement of Lemma 4 in |CHMV| . The general term of this series will be of the form (we 
adapt to our case): (a n+ i — a n )' '. In [CHMV , due to the presence of the indifferent fixed 
point, the term (a n +i — a„) decays polynomially like n~ K , say, where k > 1 depends on 
the map. In order to guarantee the aforecited summability property we have therefore to 
ask an additional assumption on k, namely n > i . We do not have such a constraint in 
our case since the length (a n+ i — a n ) decays exponentially fast. 
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(Hi) \DT(a p -i)\DT(a p _ 2 ) ■ ■ ■ DT(a' l )DT(a' ) > a", Vp < p* 
for 1 < a" < d (1)0 ) A a'. 



1 + 



\og(a"a 
log a' 



Then the first return time maps Tj and Tj n , n > 1 have the derivative 
uniformly bounded below away from a" . 

Proof. We give the proof for I and we generalize after to all the I n . We 
represent with an arrow " — > 



the evolution under T of a subset Z p C 



J, p > 1, given in (56) and (57). Consequently (bi,a ) — > (xo,a ) and 
(a' Q , b[) — > (xo, ao). In the latter case the derivative of the map coincides with 
that of T and is larger than 1, since 7\ has derivative larger than a' > 1. 
The former case follows by condition (i). For p > 1 we have: 



(b 2 , b x ) -> (a , ai) (a' , x ) , p 
->K,fl'o)-^K^o) P>3 



p-2' a p-3) ~ ^ ( a p-3> a p-4) 



(5f 



and 



(fi^&i) (ao, oi) -> K^o 
( & p-i> & p) -> K- 2 ,a. 
-^K,fl'o)->K^o) P>3 



p = 2 



(59) 



In order to get that the derivative of T p is larger than one we need: 
• in the first case 

\DT{b p ^)DT{a p ^)\DT{a' p _ 2 ) ■ ■ ■ DT(a' 2 )DT(a' 1 ) > 1 ; 



in the second case 

DTib'^lDTia^lDTia'^) 



■ ■ DT(a' 2 )DT(a' l ) > 1 



Let us suppose now that we have, for p > 1, 

\DT(ap-i)\ DT(a' p _ 2 ) ■ ■ ■ DT{a\)DT(a' Q ) > a" . 



(60) 



(61) 



(62) 



If this condition holds, then the inequality (61) follows too with the same 



uniform (in p) bound a", since, by monotonicity of the first derivative, 



DT(b ' j) > DT(a' ). In order to satisfy the inequality (60) with a lower 



bound given by a", by assuming again (62), it will be sufficient to show that 
\DT(bp-x)\ > DT(a' ), which, by monotonicity, is implied by \DT(b\)\ > 
DT(a' Q ) and this follows from assumption (ii). Therefore, we are left with 
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the proof of the validity of condition (62). By condition (iii) this holds true 



for p < p*. Moreover, since the points a' p _ 2 , a' p _ 3 , . . . , a' lies on the left of 
Xq, all the (p — 1) derivatives in the block DT(a' p _ 2 ) ■ ■ ■ DT(a' 2 )DT(a' )) are 
larger a'. On the other hand, the derivative in a p _i is surely larger than a. 



Hence, (62) holds for all p such that a (a') p_1 > a". 
Now we return to the rectangles in I n . Let us first call complete path the 
graphs given above starting respectively from (a p -i,a p ), p > 1, and ending 
in (a' ,x ) and starting from (a' ,a'_ J, p > 1 and ending in (a' ,x ). It is 
easy to check, by looking at the grammai]^] given by the arrows, that any 
subset of the rectangles in I n with first return time q > n will contain points 
whose trajectory follows a complete path, or spends some time in (xq, do). In 



any case, and by the condition (62) whose validity has been checked above, 



the derivative DT g will be strictly larger than a". ■ 

Remark 2 The assumption (i)-(iii) in the previous Lemma are easily ver- 
ified for the map investigated in t r PM$. In particular, with the values a = 
0.4603 and a' = 1.113 associated to the map and with a" ~ 1.01, the inequal- 



ity (iii) is verified for p > 9; hence we have only to check (62) for 1 < p < 8 



and this has been done, and confirmed, by a direct easy numerical inspection. 

As a consequence of the preceding results, we could apply, as in |CHMVj , 
the L.-S. Young tower theory and conclude the following statements: 

• On the induced set /, the tail of the Lebesgue measure of the set of 
points with first return bigger than n, to be more precise the quantity 
^2k>n m i x e I i T i( x ) — k}> where ti(x) denotes the first return of 



the point x into /, decays exponentially fast with n. By using (56) and 
the asymptotic values for the b n and b' n given below it is immediate 

_ n 

to find that the previous rate of decay is 0{{a') This implies the 
existence on the Borel a-algebra B ([0, 1]) of an absolutely continuous 
invariant measure [i with exponential decay of correlations for Holder 



2 Let us give the coding for the map T with the grammar which we invoked above. 
To use a coherent notation we will redefine a_o = a^; a_ p = a!, p = T 1 p a' , p > 1. 
We associate with each point x € [0, l]\x, where \ = Ui> T~ l {x Q }, the unique coding 
x = (ujq, oji, . . . ,cJ ra , . . . ), ioi G Z, where (from now on n will denote a positive integer 
larger than 1), uji = n iff T l x £ (a„_i,a„); uji = —n iff T l x £ (a_„, a_(„_i)); oji = iff 
T l x £ I. The grammar is the following (the formal symbol —0 must be intended as 0) : 

LJi = n > =>■ cu i+ i = —n ; u;,; = — n =>■ uii+i = — (n — 1) 
uji = =>■ cjj-|-i = n > (any n) 
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observables evolving under T w.r.t. p (the rate of this decay will be of 
the type a~ n , where a is possibly different from a'). 

• Since the first return maps Tj and T/ n are aperiodic uniformly expand- 
ing Markov maps, they admit invariant measures pi and pi n which turn 
out to be equivalent to Lebesgue on / and I n with densities bounded 
away from and oo |CHMVj and also Lipschitz continuous on the im- 
ages of the rectangles of the their associated Markov partition [ADj^J 
In the sequel we will show that such densities coincide, but a constant, 
with the restriction, on the inducing sets, of the density p of the invari- 
ant measure p for the map T. Now, the images of the rectangles of the 
Markov partitions are the (disjoint) sets (a' ,x ) and (xo,a ), when we 
induce over /, and the whole intervals (a' n , a! n _^) and (a n , a„+i), n > 1, 
when we induce over the rectangles in I n . Therefore we could con- 
clude that the density of the invariant measure p is a piecewise Lip- 
schitz continuous functions with possible discontinuities at the points 
a p , tip, p > 1, a ,a' and Xq. 

We now improve this last result by showing that the density is Lipschitz 
continuous over the unit interval but in the cusp point xq. We stress that 
this result will improve as well Proposition 13 in [CHMVJ. 

Proposition 3 The density p of the invariant measure p is Lipschitz con- 
tinuous and bounded over the intervals [0, 1]. Moreover, 



Proof. We work on the induced set /. The invariant measure pi for the 
induced map T/ is related to the invariant measure p over the whole interval 
thanks to the well-known formula due to Pianigiani: 



3 It is argued in |AD| that if a is a Markov partition of the standard probability metric 
space (X, B, m, T) with distance d, then Ta C <r(ct), where u{a) denotes the sigma-algebra 
generated by the partition a, and therefore it exists a (possibly countable) partition /3 
coarser than a such that a(Ta) = er(/3). Moreover, if the system is Gibbs-Markov, as in 
our case, then the space Lip^,^ of functions / : X — > K, / G :— (X) , which 
are Lipschitz continuous on each Z £ /3, is a Banach space with the norm: H/H^ ^ = 

+Dpf, where Dpf = sup Zg/3 sup^^ IJ( ^ ( ) ~^ ) fa)l ■ The space Lip^^ is compactly 
injected into L^, which gives the desired conclusions on the smoothness of the density 
as a consequence of the Lasota-Yorke inequality. Notice that in our case m in just the 
Lebesgue measure. We denote by B(I) the Banach space Lip^^ defined on /. 




(63) 




(64) 



i j=0 
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where B is any Borel set in [0, 1] and the first sum runs over the cylinders Zi 
with prescribed first return time Tj and whose union gives /. The normalizing 
constant C r = p(I) satisfies 1 = C r ^2 i Tip,i(Z i ). 

This immediately implies that by calling p the density of pi we have that 
p(x) = C r p(x) for m-almost every x £ I and therefore p can be extended to 
a Lipschitz continuous function on J as p. A straightforward application of 
formula Q gives jCHMVj : 

p,(a n -i,a n ) = C r fii(Z n+1 ) (65) 

oo 

M°n> a 'n-l) = C r Yl ^A Z P ) ( 66 ) 
p=n+2 

Let us now take a measurable B C (a' n ,a' n ); the formula above immediately 
implies that 

oo 

p(B) = C r Vi(T- {p - n) B n Z p ) (67) 

p=n+2 

Passing to the densities we have 

» oo „ 

/ p (x) dx = C r / p (x) dx (68) 

Jb p=n+2 JT~(p~")Bnz p 

We now perform a change of variables by observing that the set B is pushed 
backward p — n — 2 times by means of T^ 1 , then once by means of T^ 1 
and finally it splits into two parts according to the actions of T^ 1 and T 2 _1 . 
Therefore, 

> / p (x) dx = > > / ^ ? , ^ — dy . 

(69) 

Since B is any measurable set in (a' n , a^-i); we have for m-almost every 
x G (a! n , a n _i), 

lrp—lrp—(p—n—2) 



| T^rrn-nfrr-lrp-lrp-ip-n^) 



p=n 



X 



oo w T _i T _i ( m _2) x 

This formula does not depend on the choice of the interval (a' n , a^-i) an d 
therefore it holds for x G (0, a' ). For the cylinders (a„_i, a„) we get similarly 
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that, for m-almost any x 6 (ao, 1), 

Since p is Lipschitz continuous inside I and the inverse branches of T are 
we conclude that p can be chosen as Lipschitz continuous over the 
disjoint open intervals (0,a' ) U (a ,l). It is now useful to observe that the 



right hand sides of ( 70 ) and ( 71 ) give exactly the expression of the Perron- 
Frobenius operator associated to the first return map and whenever x is 
chosen into /. By the existence of the left (resp. right) limit of p in a' (resp. 
ao) we immediatly obtain the continuity of p in such points. We use now this 
result to prove the continuity of the density in x . We remind that such a 
density is the fixed point of the Perron- Frobenius operator, so that it verifies 
the following equation, for any x G [0, 1] : 

o(x) = P{T " {X)) + piT ^ X)) (72) 
P[ ] \DT(Tr\x)) | \DT (T 2 _1 (x)) I ' 1 } 



which gives, for x = x 



p{xo) = WW\ + WWI ' (73) 



and this proves immediately the continuity in x . 



We now observe that assumptions (52)- (55) together with the facts that 



2~2(%0 = a'p-i, Ti(a' p ) = and T\b p = T 2 b p = a p _i, allow to get easily the 
following asymptotic behaviors (for p large) for the preimages of xo (again c 
will denote a constant independent of p and that could change from a formula 
to another): 

<~ j£y-> (^p ( 74 ) 

c c 
(x - b ) ~ p ; (b p - x ) ~ — — p- . (75) 

(a') 57 («')^ 

These formulas immediately imply that for x = b p (resp. x = b') in a 
neighborhood of Xq and for p large the derivative behaves like \DT(x)\ ~ 

c (a') p( ' B ~ 1 ^ (resp. |DT(x)| ~ c(a') p( -B^^). Since p is bounded away from 
zero and infinity on /, by the preceding scalings on the growth of the deriva- 
tive near xq we have that p(x) ~ ib 1 ^ 1 for x close to and 1, which means 
that p can be extended by continuity to zero on the right side of and on 
the left side of 1. ■ 

The preceding proposition suggests the following scaling for the density. 
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Proposition 4 



p(x) = c'x a + o(x a ), x -»■ 0+; a > 0, (76) 

p(x) = c"(l - x) 6 + o((l - x) 6 ); x -> 1~ 6 > 0, (77) 

a = 6=-^-l (78) 
and t/ie constant d and c" verifying 



1 



1 . (79) 



Proof. We use again formula (72). By using for T and its two inverse 
branches the asymptotic polynomial behaviors in and 1 given in (52)-(55), 
we get at the lowest order in x in the neighborhood of 0, 

/ / /\— a— 1 a i // — b— 1 b //a /on\ 

c (a ) x+ca x=cx. (80) 

Now, suppose a < b. Then (a / ) _a_1 w 1, which implies either a' = 1 or 
a = — 1 and both cases are excluded. On the contrary, if a > b, then 
oj~ b-1 x b ~ 0, implying a = 0, which is still impossible. Hence, we necessarily 
have a = b. We now take the point x in the neighborhood of 1. By explici- 
tating Tf 1 (x) and T 2 _1 (x) with respect to x in the neighborhood of xo and 
substituting into the Perron- Frobenius equation we get, at the lowest order 
in 1 — x, 

-% + -^ = (.-^, (81) 
(1 - x) B (1 - x)^~ 

from which we obtain 

(1 -x)- 3 ^ sa (1 -x) b . (82) 

We finally conclude that a = b = — 1. Substituting this common value 
into equation (80) we finally get the expression relating the constants c' and 
c". m 

The latter relation is a good check for the validity of the shape of the 
density in and 1. 

By assuming the continuity of p in xo, we could use the value of a given 
above in terms of the map parameter B* to guess a functional expression for 
p. In agreement with the previous considerations, such an expression could 
be 

p(x) = N (7, 5) e~' yx x 5 (l — x) s , (83) 
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where, if l v (z) is the modified Bessel function of the first kind, 

>y 2 ~ e 2 

N ^ ^ = rm xAU m ' ^ 84 ) 

0rT(l + <J) 

with 5 = a, c' = N (7, 5) and c" = iV (7, 5) e~\ 

Numerical computations performed on about 10 5 values for Casimir max- 
ima allowed us to estimate the parameters describing the local behavior of 
the map listed at the beginning of this section: 

• a' ~ 1.113 , a ~ 0.4603 ; 

• B' ~ 0.3095 , B ~ 0.2856 . 

Therefore, we get B* = B' and 5 ~ 2.2258. The fit of the empirical 
stationary distribution function performed with such parameters comes out 
to be in good agreement with the functional expression for the invariant 
density (83) and the estimated value for 7 is 7 ~ 4.26. 



Figure 2: Fit of the invariant density p (x) for the map T with the function 



given in (83) 



An interesting question is to locate the maximum of the density p. Numer- 
ical investigations suggest that this maximum belongs to [a' ,xo] (see fig.2) 
depending on the parameters which define the map T. 
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3.1 Return times 



In section II and in section IV B of |PMj . the periodic orbits of the system, 
due to its invariance under R, have been empirically classified by specifying 
that the initial condition belongs to the half space containing, say, the fixed 
point Ci and the number of rotations they perform around the fixed point c 2 
(cfr. [EE] figg- 2 and lib). 

In particular, labeling as S + the portion of S laying in the half space con- 
taining Ci, it can be shown by direct inspection that fig. lib in |PMj . which 
represents the map of the set of maximum values of C in itself associated to 
periodic trajectories starting from S + after k rotations around C2, is exactly 
the graph of the induced map of T in the appropriate scale. 

Therefore, the distribution of the number of times a trajectory of the 
system, starting from S + , winds around c 2 before hitting again S + or equiv- 
alently, starting from E_, winds around c\ before hitting again £_, is the 
same of that of the random variable V xo ^ (x) , x G (xq, 1) , being the return 
time on (xq, 1) starting from x under the dynamics induced by T. In terms 
of the already constructed invariant measure /i, this probability is given by 

oo 

v(t( Xo ,i)(x) >n;xe (x , 1)) = ^ n(r ix0}1) (x) = I ; x E (x , 1)) (85) 

l=n 

oo 

— M a n-2> 0>n-l) • 

l=n 

But the sum on the r.h.s. can be computed using the corresponding expres- 
sion evaluated in ^ and we finally get 

n 

K T (xo,i)(x) > n ; x E (x , 1)) w (a') B * . (86) 

We remark that the distribution of T( xo ,i) ( x ) > x ( x 0; 1) , is the /i a.s. limit 
of the empirical distribution of the points appearing in fig. 2 of |PMj . 

We also take the occasion to remark that the average time between two 
crossing of S corresponds to the gap between the filled bands of points ap- 
pearing in fig. 2 of jPMj which has been estimated to be about 0.66 . There- 
fore, the period of the smallest periodic orbit of the Lorenz system is about 
2-0.66 in complete agreement with what predicted by the perturbation theory 
developed in |Luj and the more rigorous estimate given in |GTj . 



3.2 Statistical stability 

A slight change in the forcing term in the Lorenz equation will also change 
the shape of the associated map T and therefore the invariant density associ- 
ated to it, which will exist provided the perturped map still satisfies (52 55). 
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At the end of the section we will give two examples of such a perturba- 
tion of the forcing contribution to the Lorenz field, the first preserving the 
original symmetry of the Lorenz system, and the second breaking it. As al- 
ready remarked in the introduction, this last type of perturbation has been 
empirically shown in |CMPj to model the impact of anthropogenic forcing 
to climate dynamics of the northern hemisphere as well as the effect of the 
sea surface temperature on the Indian summer monsoon rainfall variability 

[Knq. 

Let us denote by T e the perturbed map. We show in this section that 
under suitable assumptions the density p e of the perturbed measure will con- 
verge to the density p of the unperturbed one in the L x m norm. This kind of 
property is know as statistical stability. A former paper by Alves and Viana 
|AVj , see also the succesive paper by Alves [AO , addressed the question of the 
statistical stability for a wide class of non-uniformly expanding maps. Their 
result is based on two assumptions: (i) the perturbed map belongs to an open 
neighborhood of the unperturbed one in the C k topology with k > 2 and (ii) 
the two maps are compared throughout their first return maps defined on 
the same subset where the first return maps are uniformly expanding, with 
bounded distortion and long branches. Moreover, the structural parameters 
of the perturbed map (especially those bounding the derivative and the dis- 
tortion) could be chosen uniformly in a C k neighborhood of the unperturbed 
map. The main result of those papers is that when the perturbed map con- 
verges to the unperturbed one in the C k topology then the density of the 
absolutely continuous invariant perturbed measure converges to the density 
of the unperturbed measure in the L x m norm. Here we prove the same result 
but allowing the perturbed map to be close to the unperturbed one in the 
C° topology only. We will make use of induction but, in order to preserve 
the Markov structure of the first return map, we will compare the perturbed 
and the unperturbed first return maps on different induction subsets. The 
difficulty will therefore arise in the comparison of the Perron-Frobenius op- 
erators, which will now be defined on different functional spaces. The proof 
we give is inspired by the recent work |BVj . but it contains the important 
improvement of changing the domains of inductions. Contrarily to [AV] we 
are not able to establish the continuity of the map T e t— > p e and this is surely 
due to the fact that we only require the maps C° close. On the other hand, 
discarding regularity allows us to cover a much wider class of examples; we 
believe in fact that our techniques could be used to prove the statistical stabil- 
ity for general classes of maps with some sort of criticalities and singluarities. 

Assumptions on the perturbed map 
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Assumption A T e is a Markov map of the unit interval which is one-to-one and onto 
on the intervals [0, x e p) and (x e>0 , 1], convex on both sides and of class 
C 1+tE on the open interval (0,x ei0 ) U (x efi , 1). 

Assumption B Let ||-|| denotes the C°-norm on the unit interval, then 

lim \\T e - TIL = . (87) 

Moreover, Vx G [0,1], x ^ xo, we can find e(x) such that, 
Ve < e(x), DT e exists and is finite and we have 

lim DTJx) = DT(x) . (88) 

e->0 



Furthermore, 



, DTJx) , , DTJx) , . 

lim lim ^ = lim lim ^ = 1 . 89 

x .+ x + e^o DT(x) x ^ x - e^o DT{x) K ' 



Assumption C Let us denote by Ch, e and i e respectively the Holder constant and the 
Holder exponent for the derivative of T e on the open interval (0, x efi ) fl 
(x e fl, 1); namely: \DT e (x) — DT e (y)\ < Ch, e \x — y\ w - for any x,y either 
in (0, x £j o) or in (x e>0 ,l). We assume Ch >€ and i t to converge to the 
corresponding quantities for T in the limit e — > 0. 

Assumption D Let us set d( e ,i,o) : = m f(fe e ,i,a e . ) \DT e (x)\. We assume d( e ,i,o) > 1 and that 
there exists a constant d c and e c = e (c? c ) such that, 
Ve < e c , |d(i,o) — d( e ,i,o) I < ^c- 

Remark on the notation. To simplify the notations we will set: 
w 'n = W w a' n -i)\W n = (a n ,a n+1 ) and we will denote by W' en = « in X,n-i) 
and W t>n = (a e>n , a e>n+ i) the corresponding perturbed intervals, where a' t n 
and a €) „, n > 1, are the preimages of the maximum point x £i o- We also set 

Z e>1 = U Z^ ; Z^ := « , fc' e>1 ) , Z^ := (6 e>1 , a e , ) (90) 

and 

Ze,n = Zl n U Z 2 „ ; Z\ n := (6^, b'^) , Z\ n := (b e , n , b^), (91) 

where T e b' en = T e b e ^ n = a £i „_i. The same notation will be used for the cor- 
responding unperturbed intervals. We denote by I e := (a f efi , a e ,o)\ {^e,o} the 
interval where we will induce with the first perturbed return map. From now 
on, we will denote by F the first return map of T over /, by F e the first return 
map of T e on I e and by P and P e the Perron-Frobenius operators associated 
respectively with F and F e . If t e T~ n z, where t = T^ 1 oTf^ ■ ■ -oT^z with 
%k = 1 or 2, we will call the sequence i±, . . . ,i n the signature of t relatively 
to z. 
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Remark 5 The preceding assumptions imply that the order of tangency ofT e 
in 0, x and 1 tends, in the limit e — > 0, to that ofT. This is the first require- 
ment to get again Lemma 1 for the perturbed map. The other requirement 
is expressed by Assumption D which guarantees the condition (i) in Lemma 
1. Notice that this condition cannot be deduced by assumptions (A)-(C). On 
the other hand, the assumptions (ii) and (Hi) of Lemma 1 are still valid for 
the perturbed map since we have only to control a finite number of relations 
among the corresponding derivatives. For istance, by using Assumptions B 
and C, we have 

\DT(a\) - DT e (a' e j)\ < \DT e (a[) - DT e {a'^)\ + \DT e {a[) - DT{a[)\ . (92) 

The first term on the right hand side of this inequality is controlled by the 
Holder continuity of the derivative of T t while the second one is controlled 
by the local convergence to DT of DT e in the limit e — > 0. However, we need 
some more informations for the first return maps, which are summarized in 
the following Lemma. 

Lemma 6 (i) For any n > 0, let t n and i e>n two preimages of order n of 
xo and x e $ respectively with the same signature with respect to these 
two points. Then, lim e ^ t e , n = t n . 

(ii) For any n > we have: 

lim ||T n — T n \L — . (93) 



(Hi) For any x ^ U^ =0 T k x , n > there exists e(x,n) such that, for any 
e < e(x,n), the derivative DT™(x) exists and is finite and moreover 

lim DT?(x) = DTHx) . (94) 



(iv) For any n > 1, let [u n , v n ], [u e>n , v e , n ] C [0, 1] such that u^ n ->■ u n , v €jTl ->■ 
v n in the limit e — > and T™ 1 [ Ue njVe n ],T n ] [ Un ,v n ] are injective on the re- 
spective images. Then, setting for any y G T™([u^ n , v e>n \) flT n ([-u n , v n }), 

TT {n) := CC iK^i^-T-W := (r Ik^])" 1 , (95) 

Te {n \y) ->■ T~( n \y) in the limit e ->■ 0. 

Proof. 
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(i) We prove it for n = 0, for n > 1 the proof will follow by induction. 
Suppose x t fi does not converge to x , then passing to subsequences, by 
compactness, there exists a subsequence e n and a point i / i such 
that x fn fi — > x for n — > oo. In such a point T(x) < 1 since T has only 
one maximum located at x . Now, \T tn (x ent o) —T(x)\ = |1 —T(x)\ > 0. 
We now fix a > and choose n large enough, depending on a, in such 
a way that for uniform convergence we get 

\T en (x en , ) -T{x)\ = \T en (x £n , )-T en (x)+T en (x)+T(x en , ) (96) 
-T(x en , )+T(x)\ 

<2\\T €n -T\\ + \T(x enfi )-T £n (x)\ 
<2a+\T(x €nfi )-T en (x)\ . 

In the limit n — > oo the second term on the right hand side of the 
previous inequality goes to zero by assumption B and by the continuity 
of T. We finally send a to zero getting a contradiction with the above 
strictly positive lower bound. 

(ii) The proof is standard and by induction and it uses the uniform conti- 
nuity of T n on the closed unit interval. 

(iii) We use induction again. Suppose the limit holds for n. Then we write 

\DT? +1 (x) - DT n+l (x) \ = (97) 
\DT e (T e \x))DT™(x) - DT(T n (x))DT n (x) + 
DT e {T n {x))DT?{x) - DT e {T n (x))DT r e \x)\ . 

Now, we know that: (a) x ^ Ul~lT~ k x by assumption, and also 
(b) x 7^ U£~QT~ fc a; ei o since by the induction assumption the derivative 
DT™(x) is well defined at x ^ U™ =0 T~ k x . We need to take e even 
smaller, than a certain e(x,n), to guarantee that \DT™ +1 (x)\ is well 
defined too. This is easily achieved since the preimages of xq and x £j o 
converge to each other according to signature and by choosing e small 
enough depending on x and n we could just get (a) and (b) at the same 
time and for n. We can now bound the previous expression by: 

\DT?(x)\\DT £ (T?(x)) - DT e (T n (x))\ (98) 
+ \DT e {T n {x))DT?{x) - DT(T n (x))DT n (x)\ . 

The second term converges to zero by the induction assumption. The 
first term can be bounded making use of the Holder continuity assump- 
tion on the derivative, namely 

\DT € (T?(x)) - DT e (T n (x))\ < C h>e \T?(x)) - T n (x)\^ , (99) 
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and of Assumption C assuring Ch, t and i e to converge to the corre- 
sponding quantities given for T. 

(iv) Let us set y n := T~ (n) (y) e [u n ,v n ] and y e ,„ := Tr {n \y) G [w e>n , v e ,„]. 
Suppose y et n does not converge to y n . Then, by passing again to subse- 
quences and by compactness, we can find y ^ y n such that 
lim^y^ = y. But y = T™(y ek>n ) = T-(y ek>n ) - T«(£) + T^y). 
For k going to infinity the last term tends to a value different from y 
since T is injective over [u n ,v n ], while the first difference goes to zero 
by (ii) above. 



It is clear that with the previous assumptions the map T e will admit a 
unique absolutely continuous invariant measure with density p e . This density 
will be related to the invariant density p t of the first return map F e on I e by 



the formula (70), with normalizing constant C er . Our next result will be to 



prove the statistical stability of the unperturbed density, namely 
Proposition 7 

i™ i>-*lk = o- (ioo) 

Proof. The proof is divided into two parts. The second part, which concerns 
the comparison of the invariant densities outside the regions of induction, will 
follow closely the proof of an analogous result given in |BV] , but in our case 
the proof will be easier since the quantities we are going to consider have an 
exponential tail contrarily the corresponding ones analysed in |BVj where the 
presence of a neutral fixed point forced those quantities to decay polynomially 
fast. The first part concerns the comparison of the invariant densities inside 
the regions of induction and this part is new. 

First part Let us suppose without restriction that the induction sets I = (a' , ao)\{xo}, 
— ( a ' t 0' a e,o)\{ x e,o} verify a' e g < a' Q , a e ,o < a o- in the following, to 
ease the notation we will simply write dx instead of dm(x) for the (nor- 
malized) Lebesgue measure on [0, 1] and, for any interval J C [0, 1] , we 
will set | J | := m (J) . We begin by bounding 

\p(x) -p e {x) \dx . (101) 

inh 

In footnote 3 we defined the Banach spaces B(I) and B(I t ), which are 
invariant respectively under the action of the Perron-Frobenius opera- 
tors P and P e . The densities p and p e belong respectively to these spaces 
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and they are Lipschitz continuous on the open intervals (a' , xo)U(xo, a ) 
and (a' t , x €t0 ) U (x tt0 , a ej0 ). In fact we have to consider the action of the 
Perron-Frobenius operators on a larger functional space namely that of 
functions of bounded variation. It is a standard result that the Perron- 
Frobenius operator associated to Gibbs-Markov maps with bounded 
distortion leaves invariant this space and moreover it satisfies a Lasota- 
Yorke inequality for the complete norm given by the sum of the total 
variation and the L x m norm, see for instance [B] for an account of these 
results. We denote by BV{T) and BV(I t ) the Banach spaces of func- 
tions of bounded variations defined respectively on the induction sets / 
and and by ||-|| B wn > ll'llw(J e ) ^ ne respective norms. We remark that 
the Lebesgue measure associated to this norms should be understood 
as normalized to the sets / and I e . Since the Perron-Frobenius opera- 
tors P and P e are quasi-compact on respectively BV(I) and BV(I t ), 
we know that, in the limit n — > oo, 

\\P n lj - p\\ BV{I) , (102) 

\\p:i h -M\ BV[h) ^o. (103) 

It will be important for what follows the convergence of the two pre- 
vious limits to be uniform with respect to e in the L^, and there- 
fore in the L^, norms. This is guaranteed by the results in |LS V] . 
in particular Lemmas 4.8 and 4.11. As a matter of fact, our first re- 
turn Gibbs-Markov maps fit the assumptions of the covering systems 
with countably many branches investigated in |LS V] . In particular, 
it can be proven that there exist two constants C an A such that 
||P n l/ — pWqq < CA n , where the constant C and the rate A have an 
explicit and C°° dependence on some parameters charaterizing the map 
and its expanding properties^} Therefore, given rj > we can choose n 
large enough, depending on r], and such that 

/ |p - Pe \dx = [ \p- P n lj + P n lj + P?l h - P?l Ie - Pe \dx 

< 2r] + I \P n l! - P"l /e | dx . (104) 
Jinh 

4 These constants can be explicitly computed using the Hilber metric approach. In 

particular C = (1 + a)D H e DHA ~ 2N ° A~ 2Ar ° , and A = (tanh . The integer iV insures 

that the hyperbolic diameter of the iterate P N ° of a certain cone of bounded variation 
functions is finite and bounded by Dh- In particular, a, A and D h are smooth functions 
of the quantities v and D entering the Lasota-Yorke inequality (see next footnote). 
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Let us introduce, for n > 2, 

Pn := P^lj ; p e , n := P^U (105) 

and finally p n := p n on / fl I e and p n := a n , on I e \(I fl J e ), where 
a n = \im x ^ a /+ p n (x). Notice that this right limit exists since p n is Lip- 
schitz continuous on (a' ,Xo) and moreover p n G BV(I e ) as proven in 
Section 2. We remark that the need of considering BV (I £ ) follows by 
the fact that p n could be discontinuous in x £i o- Let us rewrite the second 



term in (104) as 

/ \F n l I -P?l It \dx= [ \Pp n -P e p 6>n \dx (106) 
Jinh Jinh 



< / \Pp n - P e p n \ dx + / \P e p n - P e p €:n \dx . 
Jinh Jinh 

We now consider the term J InI \P e p n — P e p ttn \dx; by the positivity and 
the contraction in of the Perron-Frobenius operator, we have 



Pn - Pe,n\dx 

(107) 



/ \PeP~n ~ PePe,n\dx < / \P e p n - P e p e>n \dx < \ 

< / \f>n-pe,n\dx+ / \d n - p e>n \dx 

Ji e ni Jl e \(l e nl) 

= / \Ppn-l ~ P t pe,n-l\dx + / \d n - p^dx 

Jhni Jl e \(I e nI) 

< / \Ppn-l — Ptf>e,n-l\dx + m(I e \(I e fl J))(| \p n \ |oo + | \f>e ,n | |oo) 
Jl ( nl 

where the L^-norm should be understood in terms of the normalized 
Lebesgue measures respectively on / and I e . But each of these norms 
is bounded by the Banach norm and in particular for p n we have, by 
the Lasota-Yorke inequality, 

\\Pn\L < \\Pn\\ BV{I) < WP^lhvW < ^ UiWeVH) + » ■ i^) 

This last quantity, for all n large enough, is less than a constant C<i and 
the same argument also apply to Hpe.nJIoo • Moreover, setting C2 and C e ^ 



the constants bounding (108) in the unperturbed and perturbed case, 



for e sufficiently small, we have that the difference IC2 — C e ^\ is bounded 
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by a constant independent of By setting 



G r .= [ \Ppi-PtPi\dx (109) 
Jl e nl 

with I = 1, • • • , n and p\ := 1/, p t , 1 := lj e we have 

„ n 

/ \Pp n - P e p e , n \dx = J2 Gi + (n - l)C 2 m(I e \(I e n J)) (110) 
where m(I e \(I t n J)) = 0(e). 

In order to compute the term we have to use the explicit structure 
of the Perron-Frobenius operator. In particular we have 



\Pp n - P e p n \dx = 



E 



i> 



-\DF{Fr' x )\ j^\DF t {F-lx 



E 



Pn^x) 



dx < 



E 

i>i 



^(Ff 1 ^) Pn(^^) ^(i^-^) PniFjx) 



\DF(Ff 1 x)\ \DF e {F-?x)\ \DF t {F^x)\ \DF € {F$x)\ 

(111' 



Actually what we want to do is to compare the preimages of the per- 
turbed and of the unperturbed first return maps whose direct images 
are defined on cylinders with the same return times. This can always 
be done and in particular we will consider points x whose perturbed 
and unperturbed preimages are both defined. At this regard, it will be 
enough to erase from I e fl I the open interval with endpoints x e o, xq 



5 The constant v < 1 and D are in fact explicitly determined in terms of the map, 
we defer to [B] for the details. To compare with what stated in [B], we need to show 
that there exists a power no of the first return map F having the absolute value of its 
derivative uniformly larger than 2. In the case of interest, this follows easily from the 
proof of Lemma 1 by combining the Markov structure of F with the lower bound for the 
absolute value of its derivative which is uniformly larger than 1 and which is an explicit 
function of the parameters describing the local behavior of the map T, in particular a', a 
and d(i,o)- The quantities v and D are then functions of the lower bound of \DF n °\ and 
of the constant, which we denote by D', appearing in the Adler's condition. This last 
condition is equivalent to prove that T has bounded distortion and we defer to [CHMVl 
where the constant bounding the distortion is explicitly determined as a function of the 
parameters defining the map. In the present case, a simple inspection of the proof in 
[CHMV shows that such a constant is a multiple of g?(o,i)- Hence, a contribution to D' 
comes from d( .i), while the other one |CHMVj comes from the divergent behavior of the 
second derivative close to the fixed point. However, in our case, the Holder continuity 
assumption on the first derivative of the map and the exponential decay of the lenght of 
Zi makes this second contribution simply bounded by 1. 
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whose measure goes to zero in the limit e — > 0. We will prove that the 



sum in (111) is bounded uniformly in e in order to exchange the sum 



with the limit e — > 0. We remind that the perturbed and unperturbed 
induced first return maps are Gibbs-Markov and have bounded distor- 
tion and that Hpnlloo < C*2- Therefore, on each interval Z\ (resp. Z J ei ) 
i>l,j = 1,2, where F (resp. F €j j) is injective we have: 

— For any i > 1; j — 1, 2 and Vx, y E F(Zf), we have j^^Lt^j < D 1 
and Vx, y E F{Z{^ we have < £> 2 - 

— For e small enough, by the argument developed in the footnote (4), 
the difference \D 1 — D 2 \ is bounded by a constant independent of 



e. 



There exists y E Z{ (resp. Z 3 ti ) such that \DF (y)\ = (resp. 

\ i I 



This immediately implies that the first term in (j 1 1 1 [) is bounded by 

(112) 



y r^P^d X < c 2Dl y lz 



Similar bounds hold also for the other three terms in (111). We remind 



that the images of the Zi have length (xq — a' ) and the sum over the 
|Zj|'s gives the length of /. We can therefore take the limit e — > in 



(|lll|). Let us consider the first two terms in (j 1 1 1 1) 



E 

i>l 



\DF(Fr l x)\ \DF,(Fjx) 



E 



PniF-'x) 



\DF(F^ X ) 



\DF(F- l x) 
\DF e {F~lx) 
(113) 

We can bound this quantity making use of Lemma ^ part (iii) and 
part (iv) first and then by observing that the point F~ x does not 
coincide with x . Let us set w := F^ l x\ w e := F~fx and F = T l . 
Then, 



DF(F~ l x) 



DF e (Fr?x) 



i-1 

n 

m=0 



DT(T m w)DT t (T m w) 



DT € (T™w e )DT e (T m w) 



(114) 



We notice that \T m w — T™w e \ = 0(e), the intervals with endpoints 
T % w and T\w e do not contain x e ,o and their length tends to zero when 
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e vanishes. Therefore, 



i-l 

n 

m=0 



DT(T m w) 




exp 


DT e (T m w) 



i-l 

E 

m=0 



l^(y)l 



c, 



| T r, 



mm 1 1 

1 t Ho 



|T m w-T m u> £ |) 
(115) 

where y is a point between T % w and T\w t . Hence, by Assumption B, 
this term tends to 1 in the limit e — > C0 
Moreover, the other couple of terms in (111), 



E 

i>l 



< 



E 



(116) 



We remark that the function p n is a continuous extension of p to 
J e \(7 D I e ) and therefore we can rewrite 



\p n {Fr 1 x)-~ Pn {F7lx)\ = {p^F^x) - p n (F7?x)\ 



;ii7) 



where p„ is now defined on I U I e . This function is continuous on 

/ U / e \{x } and, by part (iv) of Lemma 

lim e ^ + {^(F^x) -p n (F~lx)\ = 0. 

To resume: for n larger than a certain n(rj), 



(104) <2 V + J2Gi + (n-l)0(e) ; 



:ii8i 



i=i 



each G/ is bounded uniformly w.r.t. e and tends to zero for e tending 
to zero. Therefore we can pass to the limit e — > and then r\ — Y 0. 

Second part According to the assumptions made at the beginning of the first part 
and without loss of generality, we will assume that all W e>n lie to the 

6 Actually y depends on e, y = y e . Setting y* :— lim e _>o y e , we get 

\DT e (y e ) - DT(y*)\ < \DT e (y e ) - DT e (y*)\ + \DT £ {y*) - DT(y*)\ . 

The first term on the r.h.s. can be bounded making use of the Holder continuity as- 
sumption on DT e , the second making use of Assumption B. 
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left of the corresponding W n . Therefore we have 

/ \p-p e \dx = / \p-p € \dx + / \p-pe 

J [0,1] Jl(M t JIC\W e ,l 

I \p — p e \dx 
Ji e nw{ 



dx 



+E 

oo 

+E 



1=1 



+ 



\p - p e \dx + 
w t nw e;l 

\p - p e \dx + 



f 

■J H 



\p — p e \dx 



w;nw' 



\p — p e \dx 



(119) 



w/\(w/nw e ' ,) 



The densities are given in terms of the corresponding densities of the 
induced subsets and of the multiplicative constants C r and C^ r . Hence, 
we should first compare the latter. Since they are surely smaller than 
1, we have 



< 



E 

i=l 



P 



dx 
m(I) 



dx 



m(I e 



(120) 



The same bound holds also choosing Zf (Z^) instead of Z\ [Z\^). The 
sum converges uniformly as a function of e since the norms of p 
and p e are bounded by C 2 and the lengths of the Z\ and Z x e i decay 
exponentially fast. We now show that passing to the limit e — > inside 
the sum this vanishes. At this regard we rewrite the previous bound as 



i=i 



dx 
P^^ + 



„ dx 
P- 



zfnzl. m(I) Jzl\{z}nzl.) m(I) 



(121) 



dx 



i=l 



6 m(I e ) 
2C 2 m(ZlAZ} i ) + C 2 



Pe 



dx 



zimml.) rn(I € 
1 



m(I) m(I e 



+ 



z}nz}. 



\P-Pe\ 



dx 
m(I e 



Each term in the last sum vanishes in the limit e — > 0, in particular the 
third term tends to zero by what stated in the first part of the proof. 
Moreover, by Lemma 1 and by the fact that the derivatives of the 
maps T and T e are strictly expanding in the neighborhood of xo, for 

x G (0, a ), we get 



EE 



-1^-1^-^-2) 



<c 3 



x 



(aa' log a') 



C 



4 • 



(122) 
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Furthermore, for x G (0,ao), 



V 1 



< 



min I -DTI 
(6 2 ,6i)u(&' 11 6y 



-i 



(123) 



Analogous bounds hold also for the perturbed map, so we can choose 
the constants C 4 , C 5 independent of e. Let us call p s (resp. p r ), the 
representations of the invariant density on (0, x ) (resp. (x , 1)) without 
the normalizing factor C r . By the previous bounds on the derivatives of 
T and the boundness of the densities on the induced spaces, it follows 
immediately that there exists a constant Cg such that the norms of 
p s and p r are bounded by C 6 . The same argument also holds for p e>s and 
p er and, since C 6 can be chosen independent of e, 1 1 p e s 1 1 oo ? llPerlloo < 
(\, 

We can now proceed to bound each term in (119). For the first one we 
get 

/ \p — p e \dx < \C r — C e>r \ / pdx + C e>r / \p — p e \dx , (124) 
Jinh Jinie ' Jinh 

which can be bounded uniformly in e by arguing as in the previous 
computations. For the second term (the third one can be bounded in 
the same way) we have 

/ \p — p € \dx < \C r — C e>r \ / pdx + C e>r / |/5 — p e>s \dx . 
Jmw c i ' Jmw e i ' JinWe i 

(125) 

The right hand side is uniformly bounded in e, in particular 

C e , r [ \p- Pe , s \dx < (C 2 + C 6 )m(J n W t>1 ) , (126) 
Jic\W €:1 

where vanishes m(I fl W ej i) in the limit e — > 0. 



We now consider the last sum in (119). Similar arguments allow to 



bound the remaining sum which is even easier to handle. We first have 



y\ / \p-p e \dx < v 

i=i J w[\(w;nw> tl ) l=1 



I C r C € r I / p s dx 

w{\(wir\W' e ;) 

(127) 



+C e%r / \p s - p e , s \dx 

<wi\(winw' tl ) 
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The sum is uniformly convergent as a function of e since 
^ ^e'j) c W{ and the length of such an interval decays ex- 
ponentially fast with rate independent of e. Finally, previous consid- 
erations imply that each term into the sum goes to zero in the limit 
e ->• 0. 

Finally we have 



E 

1=1 



\p — p e \dx < 



(128) 



E 

i=i 



\C r — C e 



p s dx + C t 



w.'nw' 



\p s - p e ,s\dx 



The preceding considerations also apply to the first sum in this formula 
proving this to vanish in the limit e — > 0. For the second sum we make 
use of the representations of p s and p e s in terms of the density on the 
induced space. Thus we have 



E E 



w;nw^ p=l+2k=h2 



» oo 

/ EE 

» oo 

/ EE 

J W("Kl p=I+2fc=l,2 



/ \Ps-Pe,a\dx< 
ft T -l T -l T -(p-l-2) r 



X) 



(129) 

^(T- k %- 2 1 T-} p - l - 2) x) 



Urp-lrp-lrp-ip-l-2) 



| DT'f>~ 1 (T^ T 2 T[ 



X 



lrri-lrp-ip-l-2) 



lrp-lrp-{jp-l-2) 



p(T k %-% 



x 



DTf (T~ k T~ 2 T~i 
n /^r-i^r-i^r-(p~ 1 ' 2 ) 

t J e\ 1 e,k 1 e,2 1 e,l 



X 



lrp-ip-l-2) 



P,(T~ k l T^T~t l ^x) 



\DTP- l (T k 1 T 2 1 T[ (p ^ 2) 



x 



x) 



( J £\ 1 e,k 1 e,2 1 e,l 



lrp-(p-l-2) \ 



\DTP- l (T~ 1 T 2 1 T- (p - l ' 2) x) 

Qi,i + Q21 ■ 

We further decompose Q±j as 

1 



\ DT r l {T-fc 2 xt l - 2) x)\ 



Qi,i = 



E E 



-l (rp-lrp-lrp-(p-l-2) 



inw Up t£ 2k t?, 2 \DTv-\T-^T[ 



x 



p{T k ^T^- l - 2) x) 



(130) 



-lrp-ip-l-2) 



p(T~ k X2Xi x ) + WW M - MT-^T-i 

Changing variables, setting y k := T k ~ l T 2 ~ 1 T^ p ~ l ~' 2 '- ] x and 
Ve,k ■= VeiVk) = T~ k l T~ 2 T^ p ^ 2] \T p y k ), since y k and y e>k belong to 



x) 



dx < 



dx+ 



dx 
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Z k p U Z* p , we get 

Qi,i 

But, 



E E 

k=l,2p=l+2 



P 



\p{Vk) - /5(y e ,fc) + KVe,k) - Pe{ye,k)\ d Vk ■ (131) 



oo oo 



2=1 p=2+2 



(132) 



;=i 



which is clearly convergent because the measure of Z p is exponentially 
decreasing. Moreover, by what has been shown in the first part of the 
proof, 

\p(Ve,k) - pe(ye,k)\dy k = . (133) 



lim 

e-S>0 



yk 



On the other hand, we first take e small enough to get y e ^ on the same 
side of Xq as yk and then we use the Lipschitz continuity property of 
p to conclude, by observing that y t) k tends to yk when e tends to zero, 
that also 

1 \p(Vk) ~ P(ye,k)\dy k = . (134) 



lim 

£-►0 



yk 



We now consider Q 2 ,i and show that it is uniformly bounded in e. As a 
matter of fact, 



E E 



(p-l-2) 



X 



WTp-^T^T^TI 



lrp-ip-l-2) 



X, 



X 



(135) 



1 - 



irr-(p-l-2) 



X 



DTi (T~ fc 1 T~ 2 1 T~ 1 < 



(p-l-2) 



X 



dx . 



We bound it by the sum of its two parts: the density will have bounded 
infinity norm; the sums in p over the inverses of the derivatives are 
bounded by a constant since the derivatives decay exponentially fast 
and the sums over I will be controlled by the measure of W[. Finally 
the same arguments that led to bound (114) apply also to the second 
factor in the previous expression proving it tends to zero in the limit 
e -> 0. 

This concludes the proof. ■ 

We end up our analysis considering two examples of the perturbed Lorenz 
system giving rise to perturbed versions of the map T of the kind discussed 
in this section. 
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Example 8 Let us consider a perturbation of the Lorenz field |5y obtained 
by adding the constant forcing field (0, 0, — e/3 (p + a)) , e > 0. The pertur- 
bation is easily seen to preserve the symmetry under the involution R of the 
unperturbed field. Arguing as in the first section, for e sufficiently small, the 
perturbed system will keep the same features of the unperturbed one, hence 
map T e is easily seen to satisfy (52-55) as well as Assumptions A-D. Here it 
follows the plot ofT e , for e = 0.5, and the plot of the fit of the invariant den- 
sity pn for the evolution under the maps T and T e , corresponding respectively 
to the choice of the Poincare surfaces T, + ,Y, e + , the last one being contructed 
as in the unperturbed case. 




Figure 3: T e (thick line), T (thin line). 
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Figure 4: Fit of the invariant density for the evolution under T (solid line) 
and under T e , e = 0.5, (dashed line). 



Example 9 We now consider the following perturbation of the Lorenz field 
|3[) realized by adding the field (e cos 9, e sin 9, 0) where e > and 9 E [0, 27r). 
The perturbed system is not R-invariant anymore, anyway, for e sufficiently 
small, the system will still have a saddle fixed point c e and two unstable fixed 
points c\,c\. Hence, for any 9 e [0, 2-k), we have two different T e , namely 
T+,T~, both satisfying (52 55) as well as Assumptions A-D corresponding 
respectively to the choice of the Poincare surfaces SI, which can be con- 
structed as in the unperturbed case. To obtain meaningful plots of the devia- 
tion from T of the perturbed maps as well as of the deviation of the associated 
invatiant densities from the unperturbed one, e has been set equal to 2.5 and 
9 to 70° as in KJMFj . 
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Figure 5: T (thin line), T e + (thick line to the left of T), T e (thick line to the 
right of T). 




Figure 6: Fit of the invariant densities for the evolution under T (solid line) 
and under T e + (dashed line). 
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Figure 7: Fit of the invariant densities for the evolution under T (solid line) 
and under T~ (dashed line). 
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